Photon-Number Squeezing in Circuit Quantum Electrodynamics 
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A superconducting single-electron transistor (SSET) coupled to an anharmonic oscillator, e.g., a 
Josephson junction-L-C circuit, can drive the latter to a nonequilibrium photon number state. By 
biasing the SSET in a regime where the current is carried by a combination of inelastic quasiparticle 
tunneling and coherent Cooper-pair tunneling (Josephson quasiparticle cycle), cooling of the oscil- 
lator as well as a laser like enhancement of the photon number can be achieved. Here we show, that 
the cut-off in the quasiparticle tunneling rate due to the superconducting gap, in combination with 
the anharmonicity of the oscillator, may create strongly squeezed photon number distributions. For 
low dissipation in the oscillator nearly pure Fock states can be produced. 

PACS numbers: 74.50. +r 05.45.-a 05.60.Gg 33.80.Wz 



In driven oscillator systems, depending on the type 
of driving, different nonequilibrium photon (or phonon) 
populations can be produced, and indeed for many ap- 
plications the production of specific photon number dis- 
tributions is crucial. Quantum cryptography and linear 
optical quantum computation require a reliable supply 
of single photons [TJ [2], whereas for quantum measure- 
ments it may be of advantage to use strongly squeezed 
photon distributions [3 . A well-known optical system 
where strongly squeezed states can be created is the 
micromaser [4 . In this article we describe how highly 
squeezed photon (phonon) number distributions can be 
produced in suitable superconducting quantum circuits. 

Recent experiments with such circuits, with super- 
conducting qubits coupled to electromagnetic resonators 
("circuit QED"), confirmed various concepts developed 
in the field of cavity QED but also provided important 
extensions, e.g., to the regime of strong coupling [5l El 
[71 El |9]. A Josephson qubit, ac-driven to perform Rabi 
oscillations in resonance with an oscillator, depending 
on the detuning either cools the oscillator or produces 
a laser-like enhancement of the photon numbers [TOl HI] • 
Similarly, a superconducting SET (SSET) biased at the 
Josephson quasiparticle (JQP) cycle can be used to cool 
or drive an oscillator [12 ^31 [TH [15] . Squeezing of the 
photon number distribution has been predicted, but as 
described in the literature it is only a weak effect [TTl [16] . 
Here we show that by exploiting the gap structure of the 
quasiparticle tunneling rate in combination with an an- 
harmonicity of the oscillator, strongly squeezed photon 
number states, close to a pure Fock state can be pro- 
duced (see fig. [l]). 

The system to be studied consists of a SSET coupled to 
an anharmonic superconducting or nanomechanical oscil- 
lator (see fig.[2|. The SSET is formed by a superconduct- 
ing island coupled via low-capacitance tunnel junctions 
to two superconducting leads. A gate voltage Vq shifts 
the electrostatic energy of the island and controls, to- 
gether with the transport voltage the current through 
the device. The Josephson coupling of the junctions 
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FIG. 1: A strongly squeezed distribution p{ri) of photon num- 
ber states in a Josephson j unction- L-C oscillator driven by a 
SSET obtained for the parameters ^A^g = —0.1, eV — 5.6, 
A = 2.2, g = 0.01, uj = 0.4, = 0.001, E:^ = 0.1 (all energies 
in units of Ec)^ and i^/^ — 0.004. The inset shows the energy 
dependence of the quasiparticle transition rate and energy 
differences for two transitions, A^n = E-^n+i — ^i,n — eV. 
Due to the anharmonicity of the oscillator they lie above and 
below the threshold. 



should be small compared to the charging energy scale, 
Ec = e^/2C (C is the total capacitance of the island), 
and the superconducting gap A. It leads to coherent 
Cooper pair tunneling, and even in the considered limit, 
has pronounced consequences when two charge states dif- 
fering by one Cooper-pair are nearly degenerate. In ad- 
dition, quasiparticles can tunnel incoherently (with rate 
(X 7) when the energy difference between initial and fi- 
nal states is sufficient to create the quasiparticle exci- 
tation, i.e., when it exceeds twice the superconducting 
energy gap (assumed equal for electrodes and island), 
|A£^| > 2 A. At the threshold tunneling sets in with a 
sharp step. 

The SSET is tuned to the regime of the JQP cy- 
cle [iTl [18] , where the current is carried by a combination 
of Cooper pair transfer through one junction onto the is- 
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FIG. 2: A SSET with Josephson coupling Ej and quasipar- 
ticle tunneling rate proportional to 7, coupled with strength 
g to an anharmonic oscillator. The oscillator's dissipation is 
characterized by the parameter k. To get the right sign for 
the anharmonicity we use the flux The current is deter- 
mined by the transport voltage V and the gate voltage Vg 
which is coupled to the island by the capacitance Cg- 



should be of the same order as the Josephson coupling Ej . 
The anharmonicity Q is weak, ^{v?) <C oj{n)^ where (n) 
is the average number of photons. The oscillator couples 
to the SSET with strength g. Without loss of generality 
we choose the oscillator to be at its equilibrium position 
for the island charge N = 1. 

In addition to the coherent dynamics, governed by i^o, 
the state of the system evolves due to incoherent quasi- 
particle tunneling in the SSET and due to dissipative 
processes in the oscillator. These effects will be described 
in the frame of a Liouville equation. 

We tune the gate charge close to a symmetry point, 
A^G ^ — such that only the charge states TV = 0, 1,2 
are of importance [20 . The two even states. A/" = 0,2, 
are Josephson coupled to form the basis states 

I T) = cos||iV = 0)-sin||iV = 2), 

U) = sin||iV = 0) + cos||iV = 2), (2) 

with energies 



land followed by two consecutive quasiparticle tunneling 
processes through the other junction. The energy for 
this process is provided by the voltage source. In oder to 
enhance the effect we consider in the following an asym- 
metric SSET, similar to those studied in Ref. [12 , with 
one junction (where the Cooper pair is transferred dur- 
ing the JQP cycle) having a much stronger Josephson 
coupling than the other one. 

The charge on the SSET is coupled capacitively, with 
strength characterized by ^, to an anharmonic oscillator. 
This oscillator can be realized, e.g., by a circuit combin- 
ing a capacitor, an inductor, and a non-linear element 
such as a Josephson junction, as shown in fig. [2] To 
get strong squeezing we need a positive quartic term in 
the potential energy, which can be achieved by shifting 
the minima of the inductive energy and Josephson cou- 
pling relative to each other by biasing the circuit with a 
half-integer number of flux quanta. Also nanomechani- 
cal oscillators must have an anharmonicity with positive 
sign (as studied, e.g., in Ref. [19 ) for the phonon number 
squeezing mechanism described below to work. 

The coherent dynamics of the coupled SSET and os- 
cillator is described by the Hamiltonian 

i^o = ^c(A^-A^g)'-^jcos((/)l) 

+ g{N -I) {o) + a) + ujo)a + Vt {o) of . (1) 

The charging energy depends on the number of charges 
on the island, = A^l — where A^l/r are the num- 
bers of charges which have tunnelled across the left/right 
junction and the gate charge Nq = CgVg- Cooper 
pair tunneling, assumed to dominate across the left junc- 
tion, depends on the phase which is conjugate to 
A^L, [A^L, 6*"^^] = 26*"^^. The oscillator's eigenfrequency u 



ey^ = (1 + 5Nl)Ec ± \^El^mNlEl. (3) 

Here tanx = E:s/A5NgEc, and SNq = A^q - 1 is the 
deviation from the symmetry point. In addition we con- 
sider the odd state |A^ = 1) with energy ei = EcSNq. 

We assume that the system is operated in the regime 
of vacuum Rabi oscillations, where the energy difference 
of the SSET states 5e = — is close to the oscillator 
frequency cj, and the SSET and photon number states \n) 
get strongly entangled. For weak detuning duj = u — Se 
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FIG. 3: The charge states \N = 0) and \N = 2) (here shown 
exactly at degeneracy) are coupled by the Josephson coupling 
Ej and form the basis states | j) and | |). Tuned to a reso- 
nance with the oscillator they form the eigenstates of 
the coherent part of the Hamiltonian. Also the odd states 
|l)|n) are shown. Their energy is shifted by an amount of the 
order of Ec relative to the qubit states, but to fit into the 
plot this shift is not drawn to scale. 
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FIG. 4: The energy differences A^i,n^±,n+i as a function 
of photon number n for the same parameters 6uj, Q and g 
as used in fig. [l] Also shown is a cycle of transitions, which 
increases the number of photons in the oscillator. The dashed 
arrows correspond to the rate Fj^^^^^^^^ and the solid arrows 
correspond to the rate F^^^^^^. The vertical transitions are 
due to the dissipation decreasing the number of photons. 

we can approximate the eigenstates of Hq by 

|+,n) = sin|| T)|n-l)+cos|| l)\n) , 

\-,n) = cos||T)|n-l)-sin|U)|n). (4) 

The rotation angle, tan 77 = 2g^/n/ 5E{n)^ depends on 
the effective detuning, SE{n) = Su -\- rt{2n — 1), and the 
effective coupling, g = g{] |(A^ — 1)| |). In addition, 
the states \N = l)\n) with a single excess charge on the 
island are eigenstates of the system. The energies of the 
relevant states are then given by 

Ei^n = ei -^Eoscin) , (5) 

where Eosd^) = + is the energy of the anhar- 
monic oscillator. 

The Liouville equation for the reduced density matrix 
of the composite system accounts for dissipative pro- 
cesses. Assuming the transition rates to be weak com- 
pared to the energy spacings we can use the rotating 
wave approximation and neglect off-diagonal elements. 
The probability of the system to be in the state \i) is 
then determined by the simple master equation 

j 

The transition rates Ti^j = ^f^j + ^i^j account for 
quasiparticle tunneling and for transitions caused by the 
dissipation of the oscillator. Since both have been de- 
scribed earlier, e.g. in Refs. [16l [21], it is sufficient to 
quote the results. 

The quasi-particle tunneling rate, assumed to domi- 
nate in the right junction, is given by 

rZ^ = \{j\e'1'-/'\i)\'l^,{AE,^j). (7) 



The operator e^"^^!^ transfers a single charge from the is- 
land to the right lead. (For the bias and low temperatures 
considered single charges tunnel only in this direction.) 
The relevant energy difference AEi^j = Ej — Ei — eV 
also includes the work provided by the voltage source. 
The strength and energy dependence of the rate are re- 
lated to the quasiparticle tunneling current, /qp(£^) = 
7 / dE'N{E')N{E' - E)f{E')[l - f{E' - E)], which de- 
pends on Fermi functions and the densities of states, 
N{E) = Q (\E\ - A) \E\/^/E'^ - /S?. At low tempera- 
ture we need to consider only inelastic transitions where 
energy is gained, i.e., AEi^j < 0. Due to the gap in the 
DOS the rate vanishes for AEi^j > — 2A, but jumps to a 
finite value if AEi^j < — 2A. For an anharmonic oscilla- 
tor this energy difference depends on the photon number 
n. Hence at a certain value of n the threshold may be 
passed, beyond which the transitions vanish abruptly (as 
shown in Fig. 1.) This cut-off is essential for the creation 
of a strongly squeezed photon distribution. 

The dissipation in the oscillator introduces further 
transitions, 

|l-exp([^,-^,]AT)| • 

The parameter summarizes all details of the oscillator's 
dissipation at the energy scale Ej — Ei. Most impor- 
tant are transitions with energy differences of the order 
of CO. We consider low temperatures, where only relax- 
ation processes occur. 

If the transport voltage V is large compared to the gap, 
the number of photons in the oscillator increases for posi- 
tive SNq, until a balance between driving and dissipation 
- which is proportional to n - is reached. The distribution 
of photon numbers is peaked around the average value 
^av ocj/hz^ and the Fano factor F = ((n^) — (n) )/{n) is 
slightly smaller than one. In this situation one observes 
some photon number squeezing [16 , however it is weak 
and easily destroyed by temperature. For negative de- 
tuning SNq the SSET, similar as an ac-driven qubit |TT] , 
can serve to cool the oscillator. 

As the voltage is decreased certain rates can be pushed 
beyond the threshold at AEi^j = —2 A. The rates with 
the largest energy differences AEi^ri^±,n-\-i will be cut- 
off first. These are given by 

A^i,n^±,,+i = et-ei-eF+ij^(n + l) (9) 
±^y^4g^{n^l)^[SE{n^l)]^. 

The upper branch will be cut-off first. The most inter- 
esting case is shown in fig. [4j where the cycle stops com- 
pletely beyond a certain number of photons. We will now 
describe how to tune the system to reach this situation. 

Strong squeezing requires that the energy difference 
AEi^n^_^n-\-i increases with increasing photon num- 
ber n (as is the case in fig. [4]). The condition 
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FIG. 5: (a) and (c): Average photon number (n) and Fano factor F as a function of the oscillator frequency u and the gate 
charge ANq, for eV = 5.62, g = 0.01, Q = 0.0005, Ej = 0.18, nf-f = 0.02 (all energies in units of Ec)- (b) and (d): Average 
photon number (n) and Fano factor F as a function of the frequency detuning Su and eSV for AA^o = —0.1. 



^(A^i^n^-,n+i) > requires 



n > 



(10) 



This condition is independent of n, i.e., the energy dif- 
ference either decreases or increases with n for all states. 



We can see from eq. (10) that for a negative detuning, 
(^cj < 0, an anharmonicity smaller than g is sufficient. 

Squeezing can be observed below a certain threshold 
for the transport voltage. For large n there is a limiting 
value for A£^i^n^-,n+i which has to be larger than — 2A; 
otherwise there is no point where the rate gets cut-off. 
From this we get a condition for the voltage 

^ {5E{1) - ./W+¥e{W) < eSV < -f/2n, (11) 

with eSV = eV — 2 A + ei — . The right hand inequality 
guarantees that there is a cut-off, and the left hand side 
guarantees that the cycle does not stop already at zero 
photons. If the conditions given by Eqs. (10) and (11) 
are met, the rates are cut-off at 



eSV {edV -Suj^ft) 
g^ + 2e6Vn 



(12) 



In order to have significant effects of the cut-off, we also 
have to require that nav > ^cut- 

We can optimize the system by choosing a negative 
detuning, duj < 0, precisely in a way that the system is 
in resonance at the cut-off, SE (ricut) = 0. This means 



5u; = Q[l-2{e6V)yf 



(13) 



In this case we get ncut = (e^^)^/^^- 

We solved for the stationary distribution of eq. (|6| in 
the product base of the charge states N = 0,1,2 and 
many Fock states with n < 200 sufficient to guarantee 
convergence. In fig. [s] (a) and (c) we plot the average 
photon number and the Fano factor for a fixed transport 



voltage. The maxima in the photon number and the 
minima in the Fano factor correspond to co = 5e or to 
higher order resonances. The region where the conditions 
([10| and ([n]) are fulfilled lies in the center of fig [s] (d) . 
Here we find strong squeezing. 

In fig. [T] we show the probability distribution for the 
oscillator states p{n ) = T r{p\ n){n\) for parameters which 
meet eqs. ([lO|, (11) and (13). For these parameters the 
Fano factor is particularly small, F ^ 0.01, and one can 
clearly see the effect of the cut-off. The rates are cut at 
n = 7. Therefore we have a sharp maximum in p at n = 6 
and then a sudden drop. The probabilities above n = 6 
are not exactly zero, because the numerically calculated 
quantum states allow for more none-zero matrix elements 
in the transition rates ^ than the approximations given 
by eqs. (|4|. However, these additional rates are much 
smaller than r^^^^_^ and F^^^^^^ n' ^^^^ 
by the significant drop in p for n > 6. 

Conclusion. A strongly squeezed photon number dis- 
tribution can be produced in a solid state anharmonic 
oscillator coupled to a SSET. It requires an oscillator 
with frequency in the GHz range, a positive quartic term, 
and sufficiently low dissipation, such that the inequal- 
ities uo > g > tz are satisfied. For the example 
presented in fig. [s] a Q-factor of the order of 10^ is suffi- 
cient. Nanomechanical oscillators, with Q-Factors of the 
order of 10^ have been produced [22], but reaching the 
GHz range is difficult. Circuits formed of a Josephson 
junction and L-C elements can satisfy the requirements 
and have the advantage of tunable anharmonicity and 
frequency, which in turn allows selecting the squeezed 
photon number state. Furthermore, these circuits can 
be coupled to superconducting qubits, which have been 
demonstrated to allow measuring the oscillator state as 
a state-dependent frequency shift [7 . 

Acknowledgment: We thank V. Brosco, 

M. P. Blencowe, O. Astafiev, and Y. Nakamura for 
stimulating discussion. The work is part of the EU 1ST 
Project EuroSQIP. 



[1] N. Gisin, G. Ribordy, W. Tittel, and H. Zbinden, 

Rev. Mod. Phys. 74, 145 (2002). 
[2] E. Knill, R. Laflamme, and G. J. Milburn, Nature 409, 

46 (2001). 

[3] V. Giovannetti, S. Lloyd, and L. Maccone, Science 306, 
1330 (2004). 

[4] B. T. H. Varcoe, S. Brattke, M. Weidinger, and 

H. Walther, Nature 403, 743 (2000). 
[5] A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.- 

S. Huang, J. Majer, S. Kumar, S. M. Girvin, and 

R. J. Schoelkopf, Nature 431, 162 (2004). 
[6] I. Chiorescu, P. Bertet, K. Semba, Y. Nakamura, 

C. S. P. M. Harmans, and J. E. Mooij, Nature 431, 159 

(2004) . 

[7] D. I. Schuster, A. Wallraff, A. Blais, L. Frunzio, R.- 
S. Huang, J. Majer, S. M. Girvin, and R. J. Schoelkopf, 
Phys. Rev. Lett. 94, 123602 (2005). 
[8] D. L Schuster, A. A. Houck, J. A. Schreier, A. Wallraff, 
J. M. Gambetta, A. Blais, L. Frunzio, J. Majer, B. John- 
son, M. H. Devoret, S. M. Girvin, and R. J. Schoelkopf, 
Nature 445, 515 (2007). 
[9] A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and 
R. J. Schoelkopf, Phys. Rev. A 69, 062320 (2004). 

[10] E. Il'ichev, N. Oukhanski, A. Izmalkov, Th. Wagner, 
M. Grajcar, H.-G. Wagner, A. Yu. Smirnov, A. M, 
van den Brink, M. H. S. Amin, and A. M. Zagoskin, 
Phys. Rev. Lett. 91, 097906 (2003). 

[11] J. Hauss, A. Fedorov, C. Hutter, A. Shnirman, and 
G. Schon, Phys. Rev. Lett. 100, 037003 (2008). 

[12] D. Astafiev, K. Inomata, A. O. Niskanen, T. Yamamoto, 
Yu. A. Pashkin, Y. Nakamura, and J. S. Tsai, Nature 
449, 588 (2007). 

[13] D. A. Rodrigues, and A. D. Armour, New J. Phys. 7, 251 

(2005) . 

[14] M. P. Blencowe, J. Imbers, and A. D. Armour, New 

J. Phys. 7, 236 (2005). 
[15] A. A. Clerk, and S. Bennet, New J. Phys, 7 238 (2005). 
[16] D. A. Rodrigues, J. Imbers, and A. D. Armour, 

Phys. Rev. Lett. 98, 067204 (2007). 
[17] T. A. Fulton et al, Phys. Rev. Lett. 63, 1307 (1989). 
[18] A. Maassen van den Brink, G. Schon, and L.J. Geerligs, 

Phys. Rev. Lett. 67, 3030 (1991); A. Maassen van den 

Brink, A. A. Odintsov, P.A. Bobbert, and G. Schon, 

Z. Physik B - Condensed Matter 85, 459 (1991). 
[19] R. Almog, S. Zaitsev, O. Shtempluck, and E. Bucks, 

Phys. Rev. Lett. 98, 078103 (2007). 
[20] Y. Makhhn, G. Schon, and A. Shnirman, 

Rev. Mod. Phys. 73, 357 (2001). 
[21] M.-S. Choi, F. Plastina, and R. Fazio, Phys. Rev. B. 67, 

045105 (2003). 

[22] A. Naik, O. Buu, M. D. LaHaye, A. D. Armour, 
A. A. Clerk, M. P. Blencowe, and K. C. Schwab, Nature 
443, 193 (2006). 



